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We present a complete study of the Multiparticle Biased Diffusion-Limited Aggregation 
(MBDLA) model supplemented with surface difussion (SD), focusing on the relevance and effects of 
the latter transport mechanism. By comparing different algorithms, we show that MBDLA+SD is 
a very good qualitative model for electrodeposition in practically all the range of current intensities 
provided one introduces SD in the model in the proper fashion: We have found that the correct 
procedure involves simultaneous bulk diffusion and SD, introducing a time scale arising from the 
ratio of the rates of both processes. We discuss in detail the different morphologies obtained and 
compare them to the available experimental data with very satisfactory results. We also characterize 
the aggregates thus obtained by means of the dynamic scaling exponents of the interface height, 
allowing us to distinguish several regimes in the mentioned interface growth. Our asymptotic scaling 
exponents are again in good agreement with recent experiments. We conclude by discussing a global 
picture of the influence and consequences of SD in electrodeposition. 
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I. INTRODUCTION 

Quasi-two-dimensional (quasi-2D) electrochemical de- 
position (ECD) has become one of the most widely 
studied pattern forming processes since its recognition 
as a paradigm of non-local, non- equilibrium growth pro- 
cesses [0J§). Within this general context, a great deal 
of work has been devoted in the past fifteen years to 
experimental and theoretical studies of quasi-2D ECD. 
A first group of works deals mainly with pattern for- 
mation, its main results concerning "phase diagrams" of 
morphologies [ ^0|JTT1 , ECD as Laplacian growth process 
[p~2|— |l6[| , dynamic morphological transitions fll7| , |l8| |, etc. 
All these studies aim to understanding the principles un- 
derlying the rich variety of morphologies observed, rang- 
ing from dendritic to fractal. In addition to this line of 
research, there is a second one Jl9|-p4| whose main inter- 
est is the existence of universality and scale invariance in 
the roughness of the deposits produced P],0]- From all 
this and related research, it is now believed that complex 
structures with different morphologies arise from quasi- 
2D ECD due to the interplay of different transport mech- 
anisms, such as cation diffusion, electromigration, fluid 
convection, and surface diffusion (SD) ||-|y]. However, 
the combined effect of all these factors leads to a very 
complex process, and it is becoming increasingly appar- 
ent that ECD is not well understood yet. In particular, 
the detailed role of surface diffusion (SD) is still an open 
question that hinders our understanding of both the mor- 
phologies and the scaling of ECD aggregates. 

Much of the work mentioned in the above paragraph 
has been motivated by the quest to find a universal model 
to help understand ECD phenomena. The first model for- 



mulated with that purpose was the famous computer al- 
gorithm known as Diffusion Limited Aggregation (DLA) 
p"2[ , in which a particle diffuses on a lattice and attaches 
to the growing aggregate at the place where it first hits. It 
is not difficult to observe (see and references therein) 
that this simple model represents the zero concentration, 
quasi-static limit of ECD. Therefore, its validity as a gen- 
eral description of ECD is rather restricted because it 
does not include most of the effects involved in the pro- 
cess. However, DLA has played a seminal role as a source 
of inspiration both for continuum approaches J|j|-|27j - 
which predict some high-current properties but take into 
account neither the influence of the applied voltage nor 
the electrolyte concentration — and for more sophisti- 
cated computer models, basically modifications of DLA 
(see, e.g., [28-30 and also the paragraph below), which 



are more or less phenomenological and concentrate on 
changes in morphology, thus being unable to explain the 
underlying mechanisms yielding those patterns. 

In this paper, we report on the results of detailed nu- 
merical studies of Multiparticle Biased Diffusion Limited 
Aggregation (MBDLA) supplemented with SD. 

MBDLA is a model in the family of multiparticle DLA 
models 0|5O] -B3), in which a finite number of random 
walkers, possibly with constant concentration, is intro- 
duced instead of the single walker of DLA. Thus, the 
excluded volume interaction among the walkers leads to 
several of the effects neglected in DLA. As its main ingre- 
dient, MBDLA includes, in addition, a preferential bias 
(which had been first studied in the context of single- 
particle, DLA models by Meakin Qj) of the walkers to- 
wards the cathode to mimic the electric field: In this 
form, the model was successfully introduced in mm to 
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study the influence of the applied electric field on the 
composition of magnetic, amorphous CoP alloys grown 
by ECD at constant current. The main virtue of MBDLA 
is that it is a mesoscopic model embedded on a two- 
dimensional square lattice, but it reproduces the mean 
fraction of Co and P in two-cation species ECD, as well 
as the qualitative morphology of the product electrode- 
posits. In fact, the agreement between MBDLA and ECD 
experiments pq ] is quantitative, as the electrical current 
intensity and the experiment time can be directly related 
to the simulation parameters [ ^6| . Therefore, we are con- 
fident that MBDLA is a good starting point to study the 
relevance of SD in ECD and, specifically, its influence on 
the shape of the aggregates and their dynamic scaling. 
Scaling properties of MBDLA without SD were briefly 
reported in p3[ . 

The report of our results is organized as follows: We 
describe our model in Sec. II, where a brief introduction 
to the physics and chemistry of ECD is followed by a de- 
tailed account of the rules governing MBDLA. Section III 
reports on our numerical results, such as morphological 
patterns and roughness scaling. After physically showing 
that SD has to be included, we introduce three different 
rules for SD are carefully considered and compared to 
experiments, allowing to identify the proper way to in- 
troduce SD in the model. Finally, we conclude in Sec. 
IV with a discussion of our results which will allow us to 
suggest a reasonably approximate picture of ECD phe- 
nomena. A few technical details about one the rules for 
SD are given in an Appendix. 



efficients, /i c . a their mobilities, v the fluid velocity field, 
and E the electric field along the cell. The latter is re- 
lated to cation and anion concentration via the Poisson 
equation 

V • E = —V 2 (j) = —e(z c C — z a A)/e, (2) 

where <f> is the applied potential, ez c and the 
cation and anion electric charges, respectively, and e the 
dielectric permittivity of the fluid. Generally speaking, 
matter balance across the interface leads to an interface 
velocity proportional to the flux of cations and therefore, 
in the absence of any other limiting process, proportional 
to the current density as well. In addition, except for 
the region close to the cathode, we may assume electro- 
neutrality § , which in turn implies that the cation mean 
velocity is constant. The incompressible Navier-Stokes 
equations determine the velocity v of the solvent. Fluid 
convection is always present in ECD experiments, but in 
many instances it can be small enough to be safely ne- 
glected, as has been shown in actual experiments |56],[37| . 

When the cations arrive at the cathode, they reduce 
irreversibly and an aggregate of neutral particles begins 
to grow. The particles on the surface aggregate are trans- 
ported all along it due to local chemical potential gradi- 
ents. The resulting particle current conserves the number 
of particles on the surface and is given approximately by 
by (see p8| for a detailed discussion): 

J s cx -V s /«, (3) 



II. THE MODEL 



A. Basic facts about ECD 



Prior to describing in detail what MBDLA is, and in 
order to motivate and to better understand the model 
rules, we will briefly summarize the basic physics and 
chemistry of ECD, by collecting the equations commonly 
accepted to govern its main features (see, e.g., [||](| for 
further details). Generally speaking, ECD experiments 
involve two species, named cations and anions, moving in 
an incompressible viscous fluid. In very many cases, ECD 
takes place in quasi-2D cells with parallel electrodes. The 
cations move towards the cathode and the anions towards 
the anode. The basic equations for the concentrations of 
both species are as follows: 
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where C and A are the cation and anion concentration 
respectively, D c a the cationic and anionic diffusion co- 



where J s is the particle current along the surface, k the 
interface local curvature at each site, and V s the gradient 
taken along the surface. Roughly speaking, SD tends to 
reduce the interface local curvature. Finally, we note that 
the mean concentration of charge carriers in the bath is 
constant as new cations are formed at the anode upon 
arrival of the anions M. 



B. Definition and rules of MBDLA 

In this section we will define MBDLA through its evo- 
lution rules, for which we take into account the physi- 
cal equations presented in the previous section. At this 
point, we do not consider SD, whose need will be jus- 
tified in the next Section, and consequently we post- 
pone the discussion of the rules to implement SD as 
well. Thus, MBDLA is a cellular automaton defined 
on a two-dimensional square lattice of horizontal dimen- 
sion L x and vertical dimension L y (with lateral periodic 
boundary conditions and reflective boundary condition 
at the top; for the conditions at the bottom, see below), 
in which a number of random walkers (cations) are ran- 
domly distributed with concentration c. The bottom of 
the lattice is chosen to be the cathode. We do not con- 
sider the anion dynamics, but we implicitly introduce 
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it by the creation of particles and by charge electro- 
neutrality |p9|. 

The initial condition evolves in time as follows: Ev- 
ery time step a walker is chosen and moved to one of 
its four neighboring sites with probabilities taken from a 
finite differences scheme of Eqs. and (|iq ) [jig] : prob- 
ability 1/(4 + p) to move either left, right or upwards, 
and probability (1 + p)/(A + p) to move down, i.e., to- 
wards the cathode. The parameter p is referred to as 
the bias; in galvanostatic conditions it can be quantita- 
tively related to the electric current density in the phys- 
ical system as shown in jL6|. Let us stress here that our 
present choice for the probabilities is different from that 
reported in jl6| and p3|| , but we have checked that the 
results hardly differ with those presented in this paper. 
The main reason for this new election is that, with the 
new rules, the bias p ranges from to oo, that is, from 
pure multiparticle DLA to ballistic deposition, whereas 
the rules in the mentioned references allow for a range in 
p from to 0.25, and the ballistic deposition limit cannot 
be reached (although p = 0.25 is rather close already, see 
|^6|,^3|). After a destination site has been chosen, the 
particle moves if that node of the lattice is empty; if not, 
we select another particle and repeat the destination se- 
lection procedure. Once the particle has been moved, if 
the new position has any nearest neighbor site belonging 
to the aggregate, the present position of the walker is 
added to the aggregate (the cathode or bottom bound- 
ary at the initial stage) with probability s (and is able 
to diffuse over the aggregate surface if that aggregate 
position has just one nearest neighbor belonging to it, 
see the following section); otherwise it stays there (and 
is able to move again) with probability I — s. We term 
s the sticking probability; it is related to the chemical 
activation energy the cation needs to stick to the aggre- 
gate. As particles are added to the aggregate, other ones 
are created at the top of the lattice keeping the mean 
cation concentration c constant, which in fact simulates 
an infinitely high system (experimentally this means that 
the distance between electrodes is much larger than their 
lateral dimension); consequently, the flux of particles is 
constant at every stage of the simulation. 

As we have already pointed out, the model parame- 
ters are related to the physical factors influencing the 
problem. Indeed, the choice of jump probabilities for the 
random walkers in the bath provides a recurrence relation 
which is a discretized version of the continuous equations 
( |la| ) and ffc]). Therefore, the drift velocity /i c E is pro- 
portional to the bias p. When a finite number of walkers 
is considered with concentration c, we must take into ac- 
count the excluded volume, so the effective diffusion co- 
efficient and the effective drift velocity in the simulations 
are proportional to 1 — c (in a mean field approach) [|d| . 
ft is important to note that when c — > 0, i.e., the bath 
is formed by one particle alone (as in DLA), the aggre- 
gate develops tall branches which grow at the expense of 
short ones due to screening effects. Therefore, in the low 
current limit a morphological instability appears that is 



not always present in ECD experiments. The finite con- 
centration and the hard core interaction among random 
walkers simulates the cation pressure on the aggregate, 
so c is an essential ingredient in the understanding of 
the formation of electrodeposits and to prevent these in- 
stabilities (of Laplacian character) from dominating the 
whole growth process. 

One important task is the definition of simulation time 
step. In JL6| , comparison with the experiments in |35| al- 
lowed to show that the physical time and the simulation 
time measured in number of Monte Carlos trials were 
simply proportional to each other. For this reason, we 
have stuck to the definition of the time step in as 
a Monte Carlo trial, i.e., the time needed for a particle 
to jump, both if the particle does jump or if it does not. 
Notwithstanding, we have tried other time steps defini- 
tions, such as the Monte Carlo step being defined as the 
mean time for every random walker to jump at least once, 
but the results are basically the same. Some authors de- 
fine the time step for solid-on-solid growth models as the 
mean time needed to complete an aggregate layer, but 
as we will show below, ECD electrodeposits do not grow 
with constant velocity, and therefore the mean interface 
height does not grow linearly with time. We thus believe 
that, in the ECD context, this time unit would be rather 
artificial and hence we have not used it. In fact, as we 
will show below, the work reported in this paper provides 
further evidence in favor of our choice (see the discussion 
of the experiments in [B4| in Sec. Ill B below). 



III. NUMERICAL RESULTS 

A. Morphologies 

We begin the summary of our results by discussing the 
morphologies generated by MBDLA with and without SD 
and the influence of the different rules for SD on them. 
In addition, we want to compare our computer gener- 
ated morphologies to the available experimental data. 
We take as a reference the comprehensive experimental 
work of Trigueros et al. jy]] , who reported a systematic 
experimental study of different growth regimes at con- 
stant applied voltage conditions. Their work gave rise to 
a diagram of morphologies divided into different regions 
in which similar morphologies were obtained as a function 
of the applied voltage and the electrolyte concentration. 
It is important to realize that, in galvanostatic condi- 
tions, there is no linear correspondence between voltage 
and electric current of ions, and therefore, comparison be- 
tween our morphologies and those reported by these au- 
thors can only be qualitative. No similar taxonomy work 
has been performed for constant current conditions. Al- 
though the diagram in JTl| is quite complex, it encloses 
a full variety of morphologies under the label compact. 
Some authors (f^j2l[l have studied electrodeposits sys- 
tems within this regime, and hereafter we will also refer 
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to them. Finally, a recent work by Schilardi et al. |24| ] 
provides exhaustive information on the asymptotic ECD 
regimes, which have not been considered anywhere else; 
hence, their research will also be compared to ours along 
the paper. 

As we have already mentioned, from the model per- 
spective we can compare the bias p with the electric 
current density |16|,p2[ , and c with the electrolyte con- 
centration, even though the two latter magnitudes are 
not exactly coincident, i.e., an electrolyte concentration 
equal to 0.1 M does not mean c = 0.1. We will see be- 
low that the results are not very sensitive to the specific 
value of c insofar it is not very small, and thus the dif- 
ference between actual and model concentrations is not 
very relevant. The sticking probability, s, and the diffu- 
sion parameters, namely I, A and r (or equivalently r^), 
cannot be directly tuned in an experiment, although it 
is reasonable to expect that changes in the experimental 
conditions will correspondingly modify these parameters. 
How much are they modified is something we will learn 
through our computer simulations. 



1. Bias vs. Sticking Probability without SD 

Figure [j] shows a diagram of morphologies obtained 
with < p < 5 and 0.01 < s < 0.5 without SD, with a 
particle concentration c = 0.05. We have included these 
results for two reasons: First, there has been no previ- 
ous report on MBDLA morphologies, except for a brief 
discussion in [fl6|| ; and, second, we need to discuss them 
in order to understand later what is the effect of SD on 
MBDLA morphologies. It is clear from Fig. [j] that in- 
creasing the bias or decreasing the sticking probability 
yields denser aggregates, the ones obtained for p = 
and s = 1 (bottom right) being multiparticle DLA-like 
as expected (compare to ]3l|j3^ ]). This phenomenon is 
related to the stabilizing effect of the parameters p and 
s, which can be theoretically demonstrated 0). Indeed, 
the higher the values of p, the larger the flux of particles 
reaching the interface in the direction perpendicular to 
the cathode. This reduces the probability for a cation to 
stick laterally to a branch and the screening effects due 
to the Laplacian field. On the other hand, the electric 
field combined with the reduction of the sticking proba- 
bility tends to fill the interface valleys. This first result, 
namely the fact that increasing the electric current leads 
to denser aggregates, is similar to the results reported 
by Trigueros et al. ||ll|] , who observed densification of the 
aggregates with increasing applied voltage. In particular, 
we can qualitatively compare the morphological changes 
obtained by varying the bias p for a fixed s = 0.5 in Fig. 
[j], with those provided by experimental voltage variations 
(see Fig. 2 in [jll]]). We conclude that high voltages (or in 
general, high density currents) yield denser aggregates. 
So the bias p is an essential ingredient in any realistic 
ECD model. 



As a second step in our study, we have monitored other 
relevant quantities which in turn can be experimentally 
measured, in order to obtain additional information aside 
from qualitative morphological comparisons. Figure ^ 
shows the local concentration of particles in the bath, 
still without SD, at equal time intervals. We have plotted 
the concentration profiles in the stationary regime, i.e., 
after the instability occurs (see below). Thus, the mean 
number of attached particles per unit time (or equiva- 
lently, the mean interface velocity) is constant. Leger et 
al. |^J^] have reported experimental evidence consistent 
with this stationary behavior (see e.g. Fig. 5 in j|]). We 
thus see that MBDLA agrees well with their findings, 
i.e., the stationary concentration ofparticles in the bulk 
obeys approximately the equation 



C(z) = c a + (c -c a )e-( z -*°W D % 



(4) 



where z is the vertical coordinate, zq is the interface mean 
position, c a is the concentration at the anode, cq the con- 
centration at the interface, u = /x c E, and D c is the bulk 
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diffusion coefficient. As shown in Fig 
provides a good fit of our data. In Fig. 
Eq. (J4j) (dashed line) to the simulation results, showing 
a good collapse of the bulk particle density outlines for 
different times. The small deviations close to the aggre- 
gate are due to the interface roughness. The ratio D/u 
is called diffusion length. In our fits, this length turns 
out to be about 15 lattice spacings, that is, about 2 or 
3 times the lateral width of the branches for the chosen 
parameters. This result provides another of check the 
physical validity of our model, as we can compare the 
length obtained from the fit with that taken from Ref. 
||. In this paper, the diffussion length is of order 0.5 
mm, about two times the typical branch lateral width 
(of order 1 mm), so we may conclude that the diffusion 
length obtained from our model is physically consistent. 

The inset in Fig. || shows the mean concentration front 
position zq vs. time, demonstrating that, in the station- 
ary state, MBDLA leads to a constant velocity of the 
advancing front as in the experiments. 



2. Physical Relevance of SD 

The previous subsection shows that MBDLA with- 
out SD successfully reproduces some ECD experiments, 
in particular, under galvanostatic conditions with not 
very small electric current density. However, within the 
MBDLA model it is impossible to understand the unex- 
pected compactification of aggregates in low voltage ex- 
periments pT|^9|j2l| l or the columnar-like growth found 
in other situations. Unfortunately, MBDLA aggregates 
are always ramified at low bias. In fllq] , a phenomeno- 
logical explanation of compactification was proposed by 
noticing that the reduction of s leads to more compact 
aggregates. Therefore, it was proposed there that p and s 
should be related by a monotonous function, the simplest 
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case being that of a linear relationship. With this proce- 
dure, reducing the bias leads to a corresponding decrease 
in the sticking probability, and hence to compact aggre- 
gates at low bias. However, this is an ad hoc assumption 
that cannot be experimentally tested, whereas its theo- 
retical justification is not very clear. Besides that, this 
approximation does not reproduce other morphologies, as 
that reported by Lopez-Salvans et al. fl§| or Kahanda et 
al. |2l) . In view of this, it became increasingly clear that 
there was some crucial ingredient missing in MBDLA, 
and the most obvious candidate was of course SD. 

At this point, it is instructive to consider carefully the 
work by Kahanda et al. According to their results, 
as the absolute value of the overpotential decreases the 
aggregate becomes denser, and it is formed by several 
columns which are thicker at the top than at the bot- 
tom. We interpret this as a hint on the relevance of SD: 
If, when a particle arrives at the top of a column, it dif- 
fuses along the aggregate interface, and if the diffusion 
length is shorter than the column perimeter, the particle 
will not reach the base of the pillar or another column, 
with the result of a characteristic inverted triangle struc- 
ture. The onset of similar triangle structures has also 
been reported by Pastor and Rubio 19 1. We thus came 
to the conclusion that it was necessary to include SD in 
MBDLA in order to shed further light on the nontrivial 
coupling of the different transport mechanisms. 



3. Implementation of SD in MBDLA 

We have implemented SD in MBDLA in three different 
ways, all of them starting when a particle in the bulk (the 
electrolytic solution) sticks to the aggregate but has just 
one neighbor. We have first tried two simple irreversible 
rules (other similar rules yield equivalent results, so we 
do not include them here for brevity) , named rules A and 
B, and an irreversible one, named rule C: 

Rule A: The newly incorporated particle jumps al- 
ways in the same direction, either left or right parallel to 
the cathode, until it reaches a site with at least 2 neigh- 
bors or completes I jumps. This rule is similar, but not 
identical, to the one studied in [p4| for ballistic deposition 
with surface diffusion. 

Rule B: In this second rule, we allow the particle to 
perform a random walk over the aggregate surface until it 
increases its coordination number, with a constant prob- 
ability A to be permanently stuck to its current position 
(this is the so called mortal random walker ||45|| ). 

The last rule is characterized by Arrhenius-like jump 
probabilities and, what is more important, by simultane- 
ous bulk diffusion and SD: 

Rule C: This rule allows several particles to diffuse 
simultaneously. When a particle arrives to a coordina- 
tion 1 site, it sticks and jumps to one of its two near- 
est neighbors on the aggregate with probability p n — 
exp[— Eq + (n — l)E a ], where Eo >a are adimensional ac- 



tivation energies, and n the coordination number of the 
target position. If the particle new position has 2 or 3 
neighbors, it attaches to the aggregate irreversibly. Oth- 
erwise, we label the particle as a SD particle, and we 
allow it to take further steps. Thus, we have two kinds 
of diffusing particles: Particles in the bulk, distributed 
homogeneously with concentration c; and particles that 
diffuse over the aggregate surface. With probability r 
we choose a bulk particle which evolves with its char- 
acteristic rules, and with probability 1 — r a particle on 
the surface which jumps to one of its nearest neighbors 
as we have just described for the first jump. This rule is 
close in spirit to the collective diffusion rules employed in 
studies of kinetic roughening in molecular beam epitaxy 
(MBE) f46| , and in particular to MBE models beyond 
the solid-on-solid approximation [ |47| , [l8| |. 

The main difference between rules A and B with re- 
spect to rule C is that the latter introduces a charac- 



teristic time scale Td 



while in the other cases 



diffusion is instantaneous; then, the diffusing particle is 
not affected by the overall particle dynamics. As we will 
show below, Rule C is the only one which actually repro- 
duces the influence of SD on the aggregates scaling and 
morphology. In this respect, it is important to advance 
that we have found that Arrhenius-like probabilities by 
itself are not enough to model SD: Variants of rule C 
with those probabilities and without the characteristic 
time, i.e., SD kept instantaneous, lead once again to re- 
sults similar to those of rule B. All the results presented 
were obtained with Eq = 3 and E a = 1. We have chosen 
these values to have jump probabilities smaller than 1, 
but other sets of parameters yield similar results which 
we omit for brevity. Finally, another interesting point 
is that the probabilities in Rule C allow to trivially in- 
troduce temperature in the model by simply identifying 
Eo,a — * E' a /kBT. 

Rule A, by definition, introduces a diffusion length I, 
but if I S> 1 the particle jumps practically always lead 
to an increment of its coordination, as may be seen in 
Fig. |], where some morphologies are shown for different 
values of I. The inverted triangle structure typical of the 
experiments by Kahanda et al. plfl is reproduced with 
this simple rule. Nevertheless, the tops of the pillars are 
unrealistically flat; another problem is that decreasing p 
does not lead yet to a compact aggregate regime. Rule 
A is therefore not appropriate. In the case of Rule B, 
the diffusion length is introduced indirectly by means of 
the attachment probability A (see Appendix for details). 
The mean diffusion length can be shown to be given by 
Id = 1/(2A 1 / 2 ). Morphologies obtained with this rule 
are plotted in Fig. ||. Once again, and in spite of the fact 
that rule B allows the particles to diffuse randomly over 
the aggregate, the columns developed during the growth 
turn out completely flat at the top, and the option of rule 
B was excluded as well. 

These pitfalls (and similar ones found by using 
Arrhenius-like jump probabilities, which we skip for 
brevity) led us to the conclusion that instantaneous SD 
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(or limited mobility rules, in the terminology of p6| ) is 
a too drastic one approximation for ECD. Taking into 
account that the main distinctive feature of MBDLA is 
its non-local character, interactions of particles diffusing 
along the surface with newly deposited particles are ex- 
pected to be relevant. Guided by these ideas, we propose 
rule C, which incorporates this coupling by introducing 
time scales for both bulk diffusion and SD. A sample 
of the aggregates generated by MBDLA with rule C is 
shown in Fig. || The difference with the other two rules 
is immediately apparent from the plot: This more realis- 
tic rule does induce the creation of pillars as we pointed 
out above, this time similar to those reported by Ka- 
handa et al. |2l[] and Pastor and Rubio |ll| which are 
rough at the top. Moreover, the compactification of the 
aggregates at low currents appears naturally, as can be 
noticed by following the sequence of aggregates appear- 
ing on the same row (same value of r): Decreasing the 
current leads initially to less dense aggregates, until fur- 
ther reduction of the current gives rise to more compact 
aggregates. Remarkably, there is no need to change the 
sticking probability by hand as in MBDLA without SD 
or with Rules A and B. This allows us to eliminate one 
model parameter, the sticking probability, which we take 
to be s = 1 from now on. 

So far, we have seen that, while simple SD rules provide 
good results in some solid-on-solid simulation models, the 
complex dynamics of Laplacian systems does not allow 
the particles to instantaneously diffuse; rather, we must 
allow several particles to interact before they become per- 
manently stuck to the aggregate. Roughly speaking, the 
flux of particles arriving at the aggregate defines a char- 
acteristic time t p (typically inversely proportional to the 
flux, i.e., to p). Once the particles have arrived to the 
aggregate, they diffuse until they reach a site with coordi- 
nation larger than one, or equivalently, until the particle 
meets another diffusing particle, thus forming a dimer on 
the interface which cannot move anymore. A large flux 
of particles arriving at the interface (large p) will increase 
the probability of formation of those dimers, and the par- 
ticles can hardly diffuse. The situation is not so simple 
when p is small. On one hand, the deposition mean time 
t p is large, but on the other hand, the particles hardly 
experience the applied electric current, so the probabil- 
ity of attachment to a column wall before getting to the 
bottom of the aggregate increases. Thus the Laplacian 
instability is amplified leading to a compact structure 
formed by columns and grooves. This kind of instability 
has been observed in low current galvanostatic experi- 
ments JliJ . The aggregate is therefore denser but if the 
diffusion time is not long enough the interface is unsta- 
ble. It is remarkable that this simple picture in terms of 
time scales allows to understand the relevance of SD in 
ECD experiments. 

A final important remark we would like to mention is 
that, when the diffusion probability r is about 0.99, we 
have observed some evidences of what could be a mor- 
phological transition (and the subsequent change in the 



branches) similar to those reported by Lopez-Salvans et 
al. fl7|| . However, as we want to concentrate in this paper 
on MBDLA with SD as a generic model for all regimes 
of ECD experiments, we postpone a more careful study 
of this possibility to future work, where we will pursue 
the appearance of this phenomenon for different model 
parameters (such as porr). 



4- Electrolyte concentration 

To conclude the analysis of MBDLA parameters, we 
show the effect of the electrolyte concentration, c. Figure 
exhibits the morphological changes in patterns with dif- 
ferent c values ranging from 0.01 to 0.1 for different bias 
p without SD. Note that when c — > the low current 
limit is exactly the DLA growth model pq ]. Therefore, 
we should keep a finite value of c in order to diminish the 
unavoidable DLA characteristic instability. The results 
contained in the figure allow us to conclude that, inso- 
far c is not very small, the morphologies obtained with 
MBDLA do not depend strongly on the concentration, 
and therefore the fact that there is no direct correspon- 
dence between physical and simulated concentrations is 
not a drawback of the model. 



B. Dynamic Scaling 

The previous subsection shows that the inspection of 
the morphologies is a valuable method to check the va- 
lidity and relevance of the model rules. Indeed, the un- 
realistically flat aggregates obtained with diffusion rules 
A and B unable them and motivate the investigation of 
the more realistic, non-instantaneous rule C for SD. How- 
ever, in order to exploit the main virtues of MBDLA with 
SD and to compare with other relevant models and ex- 
periments, we must take some quantitative criteria, for 
example, the analysis of the interface surface roughening. 
To this end, let us define some functions related to the 
height of the aggregate at spatial position x at time t, 
given by the scalar field h(x, t). We will also review their 
basic features before discussing MBDLA properties. 

The global width (or roughness) W(L, t) is nothing but 
the rms fluctuations of the height variable h(x,t) around 
its mean value hr,(t) — (l/L) J2 X h(x, t): 

W\L.t) = ±(j^[h{x,t)~h L (t)} 2 ^ (5) 

where angular brackets stand for ensemble average. 

Generally speaking, in many growth models, starting 
from h(x, 0) = the width satisfies the dynamic scaling 
hypothesis of Family- Vicsek |l9| : 

• \ / tP if t < L z , , a . 

W*t)~{ L a if tyyLZ \ (6) 



G 



The roughness exponent a, the dynamic exponent z, and 
their ratio (growth exponent) (3 = a/ z, identify the uni- 
versality class the model belongs to. 

In the study of kinetic roughening the height-height 
correlation function is frequently used I2I : 



C 2 (l,t) = ji^2fh{x + l,t)-h{x,t 



where ppj] 



C(l,t) 



if t</ 2 , 



(7) 



(8) 



where ai oc is the so called local roughness exponent. An- 
other important function related to the height variable h 
is the power spectrum: 



S(k,t) = (h(k,t)h(-k,t)), 



(9) 



where h(k,t) = L 1 ^ 2 J2 X [^( X ^) ~ hL(t)]exp(ikx). 
S(k, t) displays a behavior consistent with the scaling 
form ll 



S(k,t) = fc~ (2Q+1) s(fci 1/z ), 



where 



u 29 if u > 1, 
u 2 



S{U) - u2a+1 . f u <<; L 



(10a) 



(10b) 



The exponent 9 takes different values depending on the 
type of scaling exhibited by the model. For instance, 
for the so called intrinsic anomalous scaling |52[ | we have 
6 = a — aioc, whereas 6 = for Family- Vicsek scaling 
(including super-roughening, i.e., a > 1). Note that this 
implies a = ai oc . 

To apply these ideas to MBDLA characterization, a 
few remarks are in order. Although, in some cases, 
MBDLA develops ramified aggregates leading to multi- 
valued interfaces, i.e., interfaces with overhangs, it has 
been demonstrated Q that the interface of the active 
zone in DLA simulations (the aggregate sites with larger 
probability of arrival) corresponds to that constructed by 
taking the topmost site h(x, t) at every horizontal posi- 
tion x. This construction does not ensure that the mea- 
sured exponents are free of interpretations [54|], but the 
exponents are consistent with theoretical and experimen- 
tal data ||55fl . The reduction of the sticking probability 
s yields denser aggregates, and overhangs do not appear 
at any stage of the simulation for low s values. Besides 
that, if SD is present the aggregates are also more com- 
pact. In all these cases the function h(x, t) is identical 
to the aggregate outline and consequently the results do 
not have any interpretation problem. 

The main scaling features of MBDLA without SD were 
already reported in pq] . Therefore, here we will briefly 
summarize them to facilitate comparison with results in- 
cluding SD, and refer the reader to p3| for the details. 



Without SD, MBDLA displays three temporal regimes: 
At early times the global width W(L, t) features (5 = 0.5, 
this value being simply due to shot noise. This stage 
corresponds to times at which the lateral correlation 
length is of the order of the lattice spacing. After this 
noisy transient, short and large length scales are gov- 
erned by different dynamics because the bulk Laplacian 
field produces nonlocal effects (screening or shadowing 
among branches). Consequently, the local and the global 
roughness exponents, a,\ oc and a, are different and the 
interface is not self-affine. The growth exponent, /3, is 
larger than that of noise (f3 > 1/2) because some isolated 
branches begin to grow independently from each other, 
which can be understood as a signature of the Lapla- 
cian instability. As a consequence, the interface width 
grows rapidly as compared with the noise fluctuations. 
At later times, branches spread by lateral growth and 
impinge upon each other. Eventually, the system reaches 
an asymptotic regime characterized by the Kardar-Parisi- 
Zhang (KPZ) universality class |5(| exponents (a = 1/2, 
(3 = 1/3, z = 3/2). The KPZ equation is the paradig- 
matic growth model without SD, and it is given by the 
stochastic partial differential equation |36| 

^l = ^ 2 h+^(Vh) 2 + V (x,t), (11) 

where v and Ao are constants and rj(x,t) is a Gaussian 
white noise with: 

{t}(x,t))=0, (12a) 



(r)(x, t)r](x',t')) = 2D6(x - x')S(t ~ t'). (12b) 

As mentioned above, the definition of the interface 
function h(x, t) neglecting overhangs might cast some 
doubts [fi4| on the validity of the exponents reported in 
p3]| . To confirm our results, we have measured the ex- 
cess velocity produced by tilting the initial substrate and 
imposing helicoidal boundary conditions [||. The inset 
in Fig. H shows that this mean velocity is well fitted by 
a parabola, as expected for KPZ behavior. It is impor- 
tant to stress that identical results are obtained using the 
jump rules in p3fl . 

Interestingly, Schilardi et al. ]24| report experiments 
with large currents (equivalent to the large values of the 
bias p) in excellent agreement with our model. They 
observe the same three time regimes: An initial tran- 
sient with a behavior which could not be measured due 
to the resolution of the experimental device; a second 
transient with f3 > 1 characterized by the growth of iso- 
lated branches, and a third asymptotic regime at which 
the interface is characterized by KPZ exponents. A plot 
of the mean interface velocity vs. time is also given, show- 
ing a crossover from the unstable regime to the stable one 
in accordance with MBDLA predictions as can be seen 
in Fig. H The global width crosses over from the insta- 
bility (/3 > 0.5) to W(L,t) - t 1 / 3 at the time pointed 
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out by the arrow. Note that MBDLA cannot yield (3 
larger than one, because of its discrete growth rules. This 
would mean that the interface width grows faster than 
the interface mean height. Finally, the evolution of the 
morphology during the experiment is also the same in 
MBDLA and in the experiment, as seen by comparing 
Fig. ||, taken from pi] , and Fig. 10, obtained in our sim- 
ulations. 

We now consider the scaling behavior of MBDLA+SD 
for the different diffusion rules. As we have pointed out in 
the preceding section, large values of the diffusion length I 
(rule A) generate flat aggregates. This means that (3 — > 
as I — > oo at early times. Fig. 11 shows the lack of uni- 
versality in the growth exponent (3: It can be seen in this 
plot that (3 decreases with I as we expected. The same 
happens with rule B: As in the case of rule A, the growth 
exponent (3 depends strongly on the attachment proba- 
bility, A (rule B). As depicted in Fig. [n], the dependence 
is similar to that of model A since the diffusion length is 
proportional to A -1 . 

The scaling behavior in MBDLA with SD given by rule 
C is more complicated. We can recognize three different 
kinds of behavior, which we summarize as follows: 

0.05 < r < 0.25. The characteristic diffusion time 
is long, and particles diffuse rather fast along the surface 
(let us recall that they are picked with probability 1 — r at 
every Monte Carlo trial) without much interaction with 
particles arriving from the bulk, thus yielding compact 
aggregates, except if p < 0.05, because then the Lapla- 
cian field creates pillars and grooves. After a short tran- 
sient the global width grows slowly and, independently 
of the applied current, the roughness exponents are com- 
patible with those of the Edwards- Wilkinson universality 
class, whose defining equation is p7|: 



dh 



vV 2 h + r](x,t). 



(13) 



Figure [12] shows the global width collapse obtained by 
rescaling the simulation time. The plot not only shows 
the Edwards- Wilkinson growth exponent, but also the r 
independence of the results on a wide range of simula- 
tion parameters. Note that the collapsing time step is 
the one defined for MBDLA without SD divided by the 
characteristic diffusion time = r _1 . Figure [l3| shows 
the collapsed power spectrum using a = 1/2 and z = 2 
(and consequently (3 — 1/4) consistent with ( JToj ) with 
6 = for Edwards- Wilkinson exponents. It is impor- 
tant to note that this kind of dynamic scaling has been 
observed in two-dimensional ECD experiments |58|. Fi- 
nally, we have to mention that the restriction r > 0.05 
is only due to the extremely long computational times 
needed to study the model for such small values of r. 

0.3 < r < 0.7. For large p, the interface is compact 
and grows with constant velocity. The scaling is similar 
to that of the preceding case. When p — > 0, initially the 
interface is rough and the growth exponent (3 is in the 
range 0.35 — 0.40 (see Fig. |l4|). Some experiments have 
reported similar interfaces at early stages of growth pffl : 



Specifically, they obtained exponents consistent with the 
linear MBE growth model universality class (a = 3/2, 
(3 = 3/8 = 0.375 and z = 4) that is, their interfaces 
could be described by the equation (B9[ : 



^=-KV 4 h + r)(x,t). 



(14) 



Note that for this model a > 1, so the interfaces gener- 
ated with Eq. (jlj) are super-rough. In our case this short 
regime ceases when the mean interface height, h(t), is 
about 8 to 10 monolayers and the global width, W(L, t), 
is about 1. This is compatible with the referred exper- 
iments except that we do not observe the super-rough 
power spectrum. Actually, in our case the tail of S(k) 
presents a time shift at large wave vectors (Fig. [l5|) which 
is incompatible |60| ] with the behavior obtained for Eq. 
(|l4|). However, the basic phenomena, such as the value 
of the effective (3 and the onset of the instability, are in 
good agreement with the experiments. After this tran- 
sient, the aggregates are still compact and develop some 
grooves (see Fig. ^|). When these grooves appear, the 
growth exponent (3 rises dramatically due to the large 
slopes produced between grooves. Figure 16 summarizes 



all this by showing the variation of (3 with time. 

0.85 < r. Finally, when the diffusion time is short, 
three completely different situations are found as a func- 
tion of the current p. For very large p, cations become 
ballistic-ally driven to the aggregate and the unstable 
transient tends to dissappear (in fact, the p — > oo limit 
is the Ballistic Depostion discrete model, which is well 
known to belong to the KPZ universality class 
When p > 1 the aggregate grows as MBDLA without 
SD with similar parameters, except that in this case the 
aggregate mean density rises. That is, we succesively 
detect a noisy initial transient, the instability associated 
with the growing branches, and the KPZ asymptotic limit 
due to the lateral growth of the branches. The interfaces 
within the unstable regime (an example of which is shown 
in Fig. |l7|) are not self-affine but present intrinsic anoma- 
lous scaling |5l],|52| : Figure [l8] shows the power spectrum 
for r = 0.85, p = 4 and s = 1. Figure |l9| shows the col- 
lapse of the power spectrum and Fig. pTTthe collapse of 
the height- height correlation function, C(l,t), achieved 
in both cases for a = 1.78, a\ oc = 0.49, z = 2.51 and 
13 = 0.71. 

For intermediate p values (between 0.25 and 1, for al- 
most every r), the aggregate is formed by several com- 
pact thin branches which grow vertically and parallel to 
one another. In this case, the notions of a rough surface 
or dynamic scaling are meaningless. Finally, for small 
p some compact branches grow at the expenses of the 
others, so typically, one or two branches grow more than 
the others. As in the case of the parallel branches, it is 
meaningless to talk about interface roughening. 
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IV. DISCUSSION AND CONCLUSIONS 



Our first conclusion is that MBDLA is a simple compu- 
tational model which incorporates in a natural way some 
of the basic mechanisms involved in ECD experiments. 
The original model (I^j2^] was already known to be in 
good agreement with some experiments |55|. In this pa- 
per, we have provided much more evidence showing that 
MBDLA explains some of the morphological changes due 
mainly to the applied electric current and, what is more 
important, it predicts the recently observed KPZ scaling 
behavior in the high current limit (for which SD is not too 
relevant) |24j] and observed also at low currents pl] |. Be- 
fore this regime is reached, there is an unstable transient 
within which MBDLA interfaces present intrinsic anoma- 
lous scaling. We believe this type of scaling is due to SD 
not being able to communicate different portions of the 
interface fast enough, so that they grow independently 
from one another. This is analogous of the anomalous 
scaling occurrin g in the non-linear surface diffusion equa- 
tion studied in [pp. In our case, the different portions 
feature a value of the roughness exponent ai oc ~ 0.5, 
similarly to the interface subject to columnar disorder 
studied in §l]|j§|. 

Secondly, the main point of our paper is that, as 
we have seen, MBDLA without SD cannot explain low 
current experiments in which the characteristic dense 
branching aggregates of high current experiments are 
replaced by compact and column-like aggregates. Our 
working hypothesis was that the latter kind of patterns 
are due to the competition between the Laplacian field 
of the cations in the disolution and the SD current on 
the interface. Thus, our ECD model, which we wanted 
to improve as to explain, at least qualitatively, the com- 
plete ECD phenomenology, should incorporate a new rule 
for the diffusion of the adatoms attached to the aggre- 
gate. Hence, we have tried out some SD rules similar to 
those often used in growth models for molecular beam 
epitaxy [ ff6| . We have verified that instantaneous dif- 
fusion rules, namely rules that "freeze" the bulk parti- 
cles while the most recently attached particle finds its 
way through the surface, do not lead to correct results in 
the low current limit, and produce very unrealistic, flat- 
topped morphologies. We have then been forced to con- 
clude that the nonlocal character of MBDLA demands 
a diffusion rule which couples the overall cation dynam- 
ics: This is the rule we have named C. It introduces a 
characteristic diffusion time = which competes 
with the time scale related to the net flux of particles 
arriving to the interface (which, in fact, is proportional 
to the applied electric current density). With this SD 
rule, the morphologies at low, medium and high currents 
are compatible with those observed by Trigueros et al. 
Ji"l| for low, medium and high applied voltages respec- 
tively. This diffusion time Td cannot be controlled from 
the experimental point of view, but fortunately there are 
wide ranges of parameters over which the simulated mor- 



phologies hardly change, which means that the descrip- 
tion of the experiments provided by MBDLA with SD 
is robust and does not need an incontrollable parame- 
ter to be tuned. We have also compared MBDLA with 
SD with the experiments reported by Pastor and Rubio 
p0yi9|] , which characterize the product interfaces by the 
MBE exponents. MBDLA seems to reproduce the latter 
behavior for very short times and short length scales as 
can be seen in Figs. [l4| and [l^, but these results are not 
too significative, as they are not as accurate as we would 
need to make any strong claim, and could be due to the 
appearance of a characteristic short length scale. There 
is a difficulty in this respect, as MBDLA scaling is intrin- 

flyl support 



sically anomalous whereas the results in 1 20 



standard, super-rough scaling. With the data presently 
at hand we have to conclude that MBDLA with SD does 
not describe quantitatively all the aspects of the very 
low current regime, but the fact it does describe most of 
them and, above all, the compactification of the aggre- 
gates, makes us confident that MBDLA with SD is a very 
good general model for ECD. 

To conclude, we note that the model presented has the 
basic ingredients of ECD phenomena: Diffusion, electro- 
migration and surface diffusion, but for this reason, we 
have to pay a big price in terms of computational time. 
MBDLA without SD is a very time-consuming model, 
and the diffusion rules make the analysis and the sim- 
ulations a patience exercise. It has certainly been an 
improvement to find that SD rule C allows us to skip 
the sticking probability parameter, thus resducing the 
parameter space, but even then if averages of relevant 
quantities over large ensembles are required, a great deal 
of computational resources would be needed. Of course, 
this disadvantage can be removed by a careful repro- 
gramming of the algorithm, but that is another line of 
research. As our goal was to identify the most impor- 
tant factors involved in ECD, we do believe that despite 
the computing limitations of the model, MBDLA with 
SD is a powerful tool to repoduce some unclear features 
of this kind of growth experiments, and has helped un- 
derstand what are the most relevant transport properties 
and how they couple in different parameter regions. We 
hope that this work suggests further experiments to find 
out whether MBDLA with SD is the complete, general 
model for ECD or else if there are still regions which need 
separate modelling. 
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APPENDIX: SURFACE DIFFUSION RULE B 

In Rule B, surface diffusion starts when a particle 
has first arrived to the aggregate and has attached to 
it (with probability s). The particle jumps with equal 
probability to one of its two nearest neighboring sites on 
the aggregate until it increases its coordination number. 
The particle has an additional probability, A, of being 
permanently attached. This kind of particle is usually 
termed a mortal random walker p5[ . The random walk 
is performed between two absorbing boundaries, namely, 
a couple of sites with higher coordination (two or three). 
One could try to determine a priori the total number N 
of jumps the particle has to perform in each realization 
drawing such number from the probability for the parti- 
cle to take N steps on a flat line if it avoids to stick N — 1 
times and die at the iVth jump. This probability can be 
easily calculated to be given by 



Pjv(A) = A(1-A) 



JV-l 



(Al) 



However, the absorbing boundaries disallow this proce- 
dure. In any case, we have compared the simulation re- 
sults by allowing the particle to perform an actual mor- 
tal random walk, and to perform a simple random of N 
steps given by Eq. (A.1). Both results are hardly dif- 
ferent. Thus, we can approximately calculate from Eq. 



(Al) the mean and variance of the maximum number of 
jumps, given by: 



N =x> 



A ■ 



(A2) 
(A3) 



For a flat interface, the particle mean position would be 
but its variance would be 



a N = N 1 ' 2 /2= 1/(2A 1/2 ), 



(A4) 



which provides the characteristic diffusion length lr> = 
1/(2A 1 / 2 ). 
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FIG. 1. Morphologies obtained with MBDLA without surface diffusion for 256 x 400 systems with a cation concentration 
c = 0.05. Other parameters are as indicated in the figure. 
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FIG. 2. Concentration profiles for a 512 x 300 system with 
parameters p = 1, s = 1 and c = 0.1 without surface dif- 
fusion. Dashed lines represent the simulation data and solid 
lines the best fit of those data to Eq. (^). The height z is 
given is lattice spacings. 
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FIG. 3. Collapsed concentration profiles using the values 
of zo obtained from Eq. (^). Inset: The mean concentration 
front position, Zo, vs time. 




FIG. 4. Morphologies obtained with MBDLA with surface 
diffusion rule A and parameters s = 1 and c = 0.1. The size 
of the system is 256 x 400 pixels. 
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FIG. 5. Morphologies obtained with MBDLA with surface 
diffusion rule B and parameters s = 1 and c = 0.1. The size 
of the system is 256 x 400 pixels. 
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FIG. 6. Morphologies obtained with MBDLA with surface 
diffusion rule C and parameters s = 1 and c = 0.1. The size 
of the system is 256 x 400 pixels. 
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FIG. 7. Morphologies obtained with MBDLA without sur- 
face diffusion for parameters s = 1 and r = 1. Other param- 
eters are as indicated. 
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FIG. 8. Interface mean height vs time with parameters 
p = 4, s = 1 and c = 0.05 without surface diffusion. The 
arrow shows the end of the unstable regime. Inset: Mean ex- 
cess velocity, in arbitrary units, for the same parameters for 
different boundary tilts, m. Circles stand for simulation and 
the solid line is the best fit to a parabola. The dashed line 
represents the expected linear growth of KPZ type. 
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FIG. 9. in situ lateral micrographs showing the interface 
evolution from t = to t — 60 min for Ag ECD at j = 1 
mAcnT 2 in 5 x 10 -3 M Ag 2 S0 4 + 10 -2 M H 2 S0 4 + 0.5 M 
Na2S04. Taken from j2l|] with kind permission from the au- 
thors. 
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FIG. 10. Sequence of snapshots of the evolution of an ag- 
gregate grown with MBDLA without SD. Parameters arc: 
p = 0.75, s = 1, r = 1 (i.e., no surface diffusion). Times (in 
our units, see text) are (top to bottom): 35 x 10 6 , 17.5 x 10 6 , 
8.75 x 10 6 , 5.25 x 10 6 , 3.5 x 10 6 , and 1.75 x 10 6 . 
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FIG. 12. Global width vs t/r d for r = 0.05, 0.1, 0.15, 0.2 
and 0.25, with p — 0.5, s = 1 and c = 0.05. The dashed line 
is a guide to the eye, with slope 0.25. 




^ FIG. 13. Collapsed power spectrum with p = 2, s = 1, 

2 , r = 0.1 and c = 0.05 using the Ed wards- Wilkinson universal- 

10 10 10 10 10 ity class exponents at six equally spaced times from 8 x 10 6 




FIG. 11. Dependence of the growth exponent /3 on: (o) 
diffusion length I (rule A) and (□) attachment probability A 
(rule B). 
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FIG. 14. Global width vs. t/r d for r = 0.35, 0.40, 0.45, 
0.55 and 0.60, with p — 0.5, s = 1 and c = 0.1. Dashed line 
is a guide to the eye. 
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FIG. 15. Power spectrum of an interface with r = 0.5, 
p = 0.5, s = 1 and c = 0.1 at times 10 r , 2 x 10 7 , 3 x 10 7 , 4 x 10 7 , 
and 5 x 10 7 . Power spectra are anomalous at short scales. 
Dashed lines are guides to the eye. 



FIG. 18. Intrinsic anomalous power spectrum with p — 4, 
r = 0.85, s = 1 and c = 0.1. Lines correspond to interfaces at 
times 3 x 10 6 , 6 x 10 6 , 9 x 10 6 , 1.2 x 10 7 , 1.5 x 10 7 , 1.8 x 10 7 , 
2.1 x 10 7 , 2.4 x 10 7 , 2.7 x 10 7 , and 3 x 10 7 . 
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FIG. 16. Evolution of growth exponent j3 with time for 
compact aggregates with grooves. Solid line: p = 0.1 and 
r — 0.45; dashed line: p = 0.1 and r = 0.5; and dot-dashed 
line: p = 0.1 and r = 0.7. 
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FIG. 17. Dynamic evolution of the height h(x,t) with 
p — 4, r — 0.85, s = 1 and c = 0.1. Snapshots are taken 
at times 1.2 x 10 7 ,2.1 x 10 7 , and 3 x 10 7 . 
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FIG. 19. Collapsed power spectrum for the five later curves 
in Fig. |l| using: a = 1.78, z = 2.51, f3 = 0.71, and 
otioc = 0.49. Dashed lines show the slope values expected 
from Eq. (|l^) for those exponent values. 




FIG. 20. Collapsed correlation function for the same pa- 
rameters as in Fig. ^ for five equally spaced simulation times 
using: a = 1.78, z = 2.51, (3 = 0.71, and a ioc = 0.49. Dashed 
lines show the slope values expected from Eq. (Q) for those 
exponent values. 
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